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Abstract: We study quasinormal modes and scattering properties via calculation of the 
S'-matrix for scalar and electromagnetic fields propagating in the background of spherically 
and axially symmetric, traversable Lorentzian wormholes of a generic shape. Such worm- 
holes are described by the Morris-Thorne ansatz and its axially symmetric generalization. 
The properties of quasinormal ringing and scattering are shown to be determined by the 
behavior of the wormhole's shape function b{r) and shift factor <I>(r) near the throat. In 
particular, wormholes with the shape function such that b'{r) k, 1, have very long- 

lived quasinormal modes in the spectrum. We have proved that the axially symmetric 
traversable Lorentzian wormholes, unlike black holes and other compact rotating objects, 
do not allow for superradiance. As a by product we have shown that the 6th order WKB 
formula used for scattering problems of black or wormholes provides high accuracy and 
thus can be used for quite accurate calculations of the Hawking radiation processes around 
various black holes. 
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1. Introduction 

A wormhole is usually a kind of "shortcut" through space-time connecting distant regions of 
the same universe or two universes. General Relativity though does not forbid wormholes, 
requires presence of some exotic matter for their existence. A wide class of wormholes re- 
specting the strong equivalence principle is called Lorentzian wormholes, whose space-time, 
unlike that of Euclidean wormholes, is locally a Lorentzian manifold. A physically inter- 
esting subclass of Lorentzian wormholes consists of traversable wormholes through which 
a body can pass for a finite time. This last class of wormholes will be probed here by the 
passage of test fields in the wormhole vicinity and observing their response in the form 
of the waves scattering and quasinormal ringing. However, before discussing the scatter- 
ing processes, let us briefly review some basic features and observational opportunities of 
wormholes. 

The static spherically symmetric Lorentzian traversable wormholes can be modeled 
by a generic Morris-Thorne metric Q, which is fully determined by the two functions: 
the shape function b{r) and the redshift function <l>(r). The axially symmetric traversable 
wormholes can be described by generalization of the Morris-Thorne ansatz, suggested by 
Teo [||. In order to support such a wormhole gravitationally stable against small pertur- 
bations of the space-time one needs an exotic kind of matter which would violate the weak 
energy condition. The latter requires that for every future-pointing timelike vector field, 
the matter density observed by the corresponding observers is always non-negative. Thus 
for supporting a traversable wormhole one needs a source of matter with negative density 
such as, for example, dark energy [||. Wormholes might appear as a result of inflation of 
primordial microscopic wormholes to a macroscopic size during the inflationary phase of 
the Universe. If wormholes exist, they could appear as a bubble through which one could 
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see new stars. Massive, rotating wormholes surrounded by some luminous matter could 
have a visible "accretion disk" near them, similar to that surrounding supermassive black 
holes. Thus, gravitational lensing should give some opportunities for detecting wormholes 
m, General features of various wormholes were investigated in a great number of papers 
and we would like to refer the reader to only a recent few, where further references can be 
found |, 0, 0, |15|, [|, |17|. 

As we mentioned before, stability of wormholes is the crucial criterium for their exis- 
tence |18, 19, 20, 21]. In order to judge about a wormholes stability, one has to treat the 
complete set of the Einstein equations with some non vanishing energy momentum tensor 
representing the matter source. Thereby, stability depends on particular geometry of the 
wormhole metric as well as on the equation of state for matter. In this paper, we shall set 
up a more generic problem: implying the existence of some matter that provides stability 
of generic traversable wormholes, whose metric is given by some generic redshift <I>(r) and 
shape b{r) functions; we are interested in scattering processes of test fields in the vicinity 
of such wormholes. In other words, we shall consider the passage of test radiation, which 
does not influence the wormhole geometry, through wormholes of various shapes. Thus, 
our aim is to learn how the wormhole shape correlates with the scattering properties of 
various fields in its background. This is a kind of probing of a wormhole's geometry by 
test fields. 

The scattering phenomenon around wormholes can be descried by a number of char- 
acteristics, two of which, the quasinormal modes and the S-matrix, we shall consider here. 
The quasinormal modes are proper oscillation frequencies, which appear in a black hole's 
response to initial perturbation. The quasinormal modes of black holes do not depend upon 
the way of the mode's excitation but only on the black hole parameters and are, thereby, 
"fingerprints" of a black hole. It can be shown [22| that the response of a wormhole to initial 
perturbations also is dominated by the quasinormal modes. The quasinormal modes are 



defined mathematically in the same way for wormholes as they are for the black holes |22 



The S-matrix for such problems can be fully described by the wave's reflection/transmission 
coefficients. We shall show that both the quasinormal modes and reflection coefficients of 
traversable wormholes are determined by the geometry of a wormhole near its throat. The 
geometry far from the throat is insignificant for scattering in the same manner as it happens 
for black holes [23|. 

Rotating black holes are known to be spending their rotational energy on amplification 
of incident waves of perturbation [^5|. This phenomenon occurs also for various rotating 
compact bodies, such as, for example, conducting cylinders, and is called superradiance 
|24, 25, 2(:]. When considering rotating traversable wormholes, one could probably expect 
that the same superradiance should take place. However, we shall show here that rotating 
axially symmetric traversable wormholes do not allow for the superradiance. 

The paper is organized as follows: Sec. II gives some basic information about traversable 
wormholes. Sec. Ill is devoted to quasinormal modes of various traversable wormholes. 
Sec. IV is about correlation of the scattering coefficients with a wormhole's shape. Sec. V 
considers scattering around rotating traversable axially symmetric wormholes. Finally, in 
the Appendix we demonstrate the utility of the 6th order WKB formula for calculations 
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of the graybody factors of black holes, which gives quite accurate "automatized" formulas 
for calculations of the intensity of the Hawking radiation. 



2. Basic features of traversable wormholes 

Static spherically symmetric Lorentzian traversable wormholes of an arbitrary shape can 
be modeled by a Morris-Thorne ansatz 

^ r 

Here $(r) is the lapse function which determines the red-shift effect and tidal force of the 
wormhole space-time. The $ = wormholes do not produce tidal force. The shape of a 
wormhole is completely determined by another function b{r), called the shape function. 
When embedding the wormhole metric into the Euclidian space-time with cylindrical co- 
ordinates on it, the equation that describes the embedded surface on the equatorial slice 
e = it/2 is 

/ r \"2 

- = ± — - - 1 . (2.2) 



dr \b{r) 

The throat of a wormhole is situated at a minimal value of r, rmin = bo- Thus the coordinate 
r runs from rmin until spatial infinity r = oo, while in the proper radial distance coordinate 
dl, given by the equation 

- ± 1 - ^ , (2.3) 



dr 

there are two infinities / = ±00 at r = 00. 

From the requirement of the absence of singularities, <I>(r) must be finite everywhere, 
and the requirement of asymptotic flatness gives <I>(r) — t- as r — )■ cxd (or / — )• iLcxd). The 
shape function 6(r) must be such that 1 — b{r)/r > and b{r)/r — )• as r — )■ 00 (or 
/ — >• ±00). In the throat r = b{r) and thus 1 — b{r)/r vanishes. The metric of traversable 
wormholes does not have a singularity in the throat and the traveler can pass through the 
wormhole during the finite time. Some thermodynamic properties of traversable wormholes 
have been recently considered in p^ ]. 

As an illustrative example, we shall consider here two-parametric family of wormholes 
which are described by the following shape and redshift functions: 

b{r) = bl~'^r'', q < 1 $(r) = -L, p > or $(r) = (2.4) 

The shapes z(r) are plotted for the above shape function in Fig. 1. The larger q corresponds 
to a smaller slope of the curve z{r) respectively horizontal axe z, that is a longer spatial 
extension of the throat regime. For all wormholes z'{r) = 00 at the throat r = bo, while for 
the ultimate case g — )• 1 — 0, z'{r) grows without bound at all points r, thus extending the 
spatial region of the throat dominance. 

In principle we could choose other forms of the function b{r) than those we mentioned. 
This would not change any of the conclusions of this work, because the natural generic form 
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Figure 1: Wormhole shapes z(r); b{r) = 6q 'r'', $ = 0. 

of the wormhole shape (when embedded on the equatorial slice) is some axially symmetric 
smooth surface with monotonically increasing section when moving away from the throat. 
This generic form of a wormhole can be represented by the functions ( |2.4D or by some other 
functions satisfying the above mentioned general requirements for traversable wormholes 
without any qualitative changes of the scattering characteristics. 

3. Quasinormal modes 

The wave equations for test scalar <I> and electromagnetic fields are 

(5''^V=5<J>,m),- = 0> (3-1) 



- A„,,)<7"'^g'^'' = 0. (3.2) 

After making use of the metric coefficients (|2.1| ), the perturbation equations can be reduced 

to the wavelike form for the scalar and electromagnetic wave functions and (see for 
instance p3| ): 

^ + {co^-V,{r))^i = 0, dn = {A{r)B{r))-'/'dr, (3.3) 
where the scalar and electromagnetic effective potentials have the form: 

V; = ^(^)^ii±ll + l(^A{r)B'ir) + A'{r)B{r)), (3.4) 
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Figure 2: The effective potential for the scalar field £ = as a function of the tortoise coordinate; 
b{r) =bo,^ = 0. 
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Here A(r) = e^*'-'"-', B{r) = l — b{r)/r. The effective potentials for the wormholes described 
by the metric ( ^.1| ) , ( p.4| ) in terms of the tortoise coordinate have the form of the positive 
definite potential barriers with the peak situated at the throat of a wormhole (see Fig. 
|2[). The whole space lies between the two "infinities" connecting two universes or distant 
regions of space. According to ||2^, the quasinormal modes of wormholes are solutions of 
the wave equation (|3.3|) satisfying the following boundary conditions: 



ibcxD. 



(3.6) 



This means the pure outgoing waves at both infinities or no waves coming from either the 
left or the right infinity. This is a quite natural condition if one remembers that quasi- 
normal modes are proper oscillations of wormholes, i. e. they represent a response to the 
perturbation when the "momentary" initial perturbation does not affect the propagation 
of the response at later time. We shall write a quasinormal mode as 



(3.7) 



where ojRe is the real oscillation frequency and ujirn is proportional to the decay rate of a 
given mode. 
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It is remarkable that the effective potential as a function of the r coordinate (see 



formulas 3.5, 3.4) from 60 to infinity has a breaking point at the throat. Though, there 
is a smooth behavior at the throat respectively the r* coordinate, so that the left and 
right derivatives of the effective potential in its maximum (throat) are equal. Therefore for 
finding the quasinormal modes and the S-matrix one can use the WKB formula suggested 
by Will and Schutz [^] and developed to higher orders in |29|. The same formula was 



effectively used in a number of works for quasinormal modes of black holes [30| as well as 
for some scattering problems 

The 6th order WKB formula reads 

-^-Y^A, = N + l, iV = 0,l,2..., (3.8) 



2' 

V --VU i=-2 

where Q = uj'^ — V and the correction terms Aj were obtained in , . Here Qq means 
the i-th derivative of Q at its maximum with respect to the tortoise coordinate r^. 

The 6th order WKB quasinormal modes are shown in Figs. 3-6. In Fig. 3 one can see 
that as q increases from zero to unity, the real part of uj for the scalar field is monotonically 
decreasing approaching some constant. For the electromagnetic field w/je monotonically 
increases with q reaching some constant at g = 1. The damping rate (Figs. 3-6) is 
monotonically decreasing as q increases approaching asymptotically zero for q = 1. This is 
a remarkable result because in the limit of the cylindrical form of geometry near the throat 
g = 1 the quasinormal modes reduce to the pure real, nondecaying, quasinormal modes. 
Such nondecaying pure real modes appear for instance in the spectrum of the massive 
scalar field of Schwarzschild [p3] and Kerr black holes and are called quasiresonances 



1 34]. Let us note that as we neglected the behavior of the effective potential far from the 
throat, the approaching of the regime of pure real modes is apparently never exact and in 
practice we can talk about long-lived modes when approaching the limit b'{r = bo) = 1. 

We can reproduce the above numerical results for higher £ in the analytical form. Using 
the first order WKB formula for the general effective potential V with arbitrary b{r) and 
$(r) and then expanding the result in powers of l/£ we obtain a relatively concise formula 



bo V 2; V260 



+ ^ TM^^r^lKWM^- (3.9) 



In the particular case of no-tidal force wormholes ($ = 0), we have 

1 1\ . 1 / 1' 



UJ 



bo V 27 V 2 



+ - n+- yjl-b'ib^), (3.10) 



These formulas are valuable because they usually give good accuracy already for moderately 
low multipoles I = 2,3,4, .. .. The above formulas let us make the two main conclusions: 
first is that the quasinormal modes are indeed determined by the behavior of the shape 
function and shiftfunction only near the throat, and second, the limit b'{r = 69) = 1 gives 
pure real nondamping modes, similar to the standing waves of an oscillating string with 
fixed ends. 
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Figure 3: Real (right panel) and imaginary (left panel) parts of quasinormal modes for scalar field 
for £ = (red, bottom), I = 1 (blue), 1 = 2 (green, top); h{r) = bl^'^r'^, $ = 0. 




Although we considered here only scalar and vector quasinormal modes, a qualitatively 
similar phenomena could be expected for gravitational perturbations at least for relatively 
large values of I because the centrifugal part of effective potential, which dominates at 
larger i, is the same for fields of any spin. Though, as we mentioned earlier, gravitational 
modes can give surprises in the quasinormal spectrum due to the exotic character of matter 
supporting a wormhole. 



4. 5'-matrix 

In this section we shall consider the wave equation ( |3.3D with different boundary conditions, 
allowing for incoming waves from one of the infinities. This corresponds to the scattering 
phenomena. The scattering boundary conditions for ( |3.3| ) have the following form: 

^. = Te"*'^'^*, r, ^ -DO, ^ ■ ' 

where R and T are called the reflection and transmission coefficients. 

The above boundary conditions ( [4.1| ) are nothing but the standard scattering boundary 
conditions for finding the S-matrix. The effective potential has the distinctive form of the 
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Figure 5: Real (right panel) and imaginary (left panel) parts of quasinormal modes for the Maxwell 
field for ^ = 1; b{r) = bl~'^ri, $ = l/r^. 



1 . 5l 



-0.5 



0.5 




Figure 6: Real (right panel) and imaginary (left panel) parts of quasinormal modes for scalar field 
for e=l; 6(r) = feJ^V"?, $ = l/r^ 



potential barrier, so that the WKB approach ||2^ can be applied for finding R and T. Let 
us note, that as the wave energy (or frequency) to is real, the first order WKB values for 
R and T will be real ^ and 

|rp + |i2|^ = l. (4.2) 

For a;2 ^ Vq, we shall use the first order beyond the eikonal approximation WKB 
formula, developed by B. Schutz and C. Will (see pM]) for scattering around black holes 



+ u^c^Vo, (4.3) 

where K is given in the Appendix. 

After the reflection coefficient is calculated we can find the transmission coefficient for 
the particular multipole number i 

lAef = 1 - = \T,\\ (4.4) 
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Figure 7: Transmission coefScients for scalar (left panels) and electromagnetic (right panels) fields 
e ^ 1 (top panels) and £ = 2 (bottom panels), $ = 0, b{r) = bl^'^r'^, q = -I (violet), q = -1/2 
(cyan), q = (blue), q = 1/2 (green), q — 3/4 (yellow), q — 7/8 (red), q ~ 15/16 (magenta). As 
g — > 1 the transmission coefficient approaches the theta function 0{uP'h^ — £{£ + 1)). 

From the Figs. 7 and 8 one can see that in the limit b'(bo) = 1, the transmission 
coefficient approaches asymptotically the theta function, which is dependent on the values 
of the shift function $(6o) (see Fig. 9) and the shape function bo on the throat, 



1. 



(4.5) 



In the Appendix one can see that the application of the higher order WKB approach 
allows us to achieve high accuracy in calculations of the transmission coefficients, and, as 
a result, of the intensity of the Hawking radiation. 

5. Rotating axisymmetric traversable wormholes 

The asymptotic forms of the solutions near the black hole and at the spatial infinity are 

$(r,) = e-''^'* + Re+''^'* , r ^ oo, (5.1) 
qj(r^) =Te-'^"^'^>'', r^r+. (5.2) 

Here R is called the amplitude of the reflected wave or the reflection coefficient, and T is 
the transmission coefficient. If \R\ > 1, that is 



> 1, 



(5.3) 
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Figure 8: Transmission coefficient for scalar (left panels) and electromagnetic (right panels) 
fields £ = 1 (top panels) and ^ = 2 (bottom panels), $ = b{r) = fej"''r«, q = -1 (violet), 
q = -1/2 (cyan), q = (blue), q = 1/2 (green), q = 3/4 (yellow), q = 7/8 (red), q = 15/16 (ma- 
genta), q — 31/32 (orange). As g ^ 1 the transmission coefficient approaches the theta function 




Figure 9: Transmission coefficient for scalar (left) and electromagnetic (right) fields i — 1, b{r) = 
6o, $ = {^Y: p = -1/2 (magenta), p = -1/4 (red), p = (green), p = 1/4 (blue), p = 1/2 (cyan). 
The more sloping line corresponds to the lower value of p. 



then, the reflected wave has larger amplitude than the incident one. This amplification 
of the incident wave is called the superradiance and was first predicted by Zel'dovich 



The superradiance effect for Kerr black holes was first calculated by Starobinsky ||2^, 26 1. 



-10- 



The process of superradiant amplification occurs due to the extraction of rotational energy 
of a black hole, and therefore, it happens only for modes with positive values of azimuthal 
number m, that corresponds to "corotation" with a black hole. Some aspects of super- 
radiance of four dimensional black holes of astrophysical interest were considered in [^] . 

Let us now see if this superradiant amplification is possible for traversable wormholes. 
The line element of the stationary, axially symmetric traversable wormhole can be written 
as @ 

ds' = -e^^dt^ + + r^K\de^ + sin2 - Codt)^), (5.4) 



1 - b/r 

where b, K, and cu, being functions of r and 9, are chosen in such a way that they are 
regular on the symmetry axis 6 = 0, vr. 

When b, and O are functions of r and K = K[9) we can separate the variables in 
the equation of motion for the test scalar field pf 



Using the ansatz 



^d^ig^'^^d^cj,) = 0. (5.5) 



we find the wavelike equation 

where \i^Tn is the angular separation constant, which depends on the azimuthal number m 
and the multipole number I > \m\. In the particular case, ii K = 1, S{6) are spheroidal 
functions, and Xi^m = ^(^ + !)■ 

The tortoise coordinate is given by 

dr 

dr^ = - 



The metric describes two identical, asymptotically flat regions joined together at the 
throat r = b > 0. Therefore, the boundary conditions for the radial part of the function at 
the infinities are symmetric 

where ^}{u}, m) = uj + mu){r = oo). 

The tortoise coordinate r^, maps the two regions r > b onto (—00,00) and for the 
scattering problem we have the following boundary conditions: 

^ = g-in(<^,m)r. ^ _^giC{^,m)r, ^ ^ OO, 
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By comparing the Wronskian of the two hnearly independent solutions of the waveUke 
equation ^ and ^* at the boundaries, we find that 

|T|2 + |i?|2 = l 

for any possible effective potential. 

This relation does not depend on the particular form of the effective potential and the 
function uj and can be also found for other test fields and gravitational perturbations. Thus, 
we conclude that there is no supcrradiance for axially symmetric traversable wormholes. 
Because of the symmetry of the regions at both sides of the wormhole, the energy, extracted 
from the wormhole when a wave falls inside it, is absorbed by the wormhole on the other 
side. 

6. Discussions 

We have considered the propagation of test scalar and Maxwell fields in the background of 
the traversable spherically symmetric and axially symmetric wormholes of a generic form. 
The majority of our conclusions are independent of the concrete shape of a wormhole. 
We have found that when the shape function 6(r) approaches the regime 6'(r) = 1 at 
the wormhole's throat b = bo, the quasinormal modes asymptotically approache pure real 
nondecaying modes, meaning the appearance of long-lived modes. The latter might give 
much better observational opportunities. The reflection coefficient approaches the step 
function which depends on the values of the shape and redshift function at the throat. We 
have proved that independently on their particular forms, traversable axially symmetric 
wormholes do not allow for superradiance. We were limited here by wormholes whose shape 
is symmetric with respect to their center (throat). Spherically (but not axially) symmetric 
wormholes are also necessarily symmetric with respect to the throat. 
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Appendix: Calculation of the reflection coeflicients for generic black and 
wormholes by the 6th order WKB formula 

In order to estimate accuracy of the WKB formula for the scattering problem, we shall 
consider the well-known Schwarzschild metric 

dr"^ 1 
ds'^ = -f{r)dt^ + + r'^{de'^ + sin^ Odct)'^), f{r) = 1 - -■ (6-1) 

The scalar field equation, 

U<t) = 0, 
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Figure 10: On the left figure we show the graybody factors for the scalar field £ = I in the 
Schwarzschild background found with the help of the WKB formula (red line) and the accurate 
values found by the shooting method (blue line) . The right figure shows the difference between the 
WKB formula and the accurate results. 



dtdcj 




Cl) 



Figure 11: Energy emission rate of the Schwarzschild black hole for the scalar field found with the 
help of the WKB formula (red, bottom) and the accurate values obtained by the shooting method 
(blue, top). 

can be reduced to the wavelike equation ( |3.3D with the effective potential 
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The reflection coefficient, given by the WKB formula, is 



i?= (l + e-2^-^)-5, (6.2) 



where 



Here Vq is the maximum of the effective potential, Vq is the second derivative of the effective 
potential in its maximum with respect to the tortoise coordinate, and Aj are higher order 
WKB corrections which depend on up to 2ith order derivatives of the effective potential 



at its maximum [29|. 



We compare the graybody factors found with the help of the WKB formula with the 
graybody factor calculated accurately by fitting the numerical solution of the wavelike 
equation (see Fig. 10). In Fig. ^ we show the energy emission rate for the scalar field. 



There one can see that the accuracy of the 6th order WKB formula is exceptionally good 
for not very small values of cj giving us at hand an automatic and powerful method for 
estimation of the Hawking radiation effect for various black holes. 

The automatized WKB formula for the quasinormal modes and reflection coefficient 
in the Mathematica(§) Notebook format and examples of using the WKB formula for 
Schwarzschild black hole can be downloaded from http://fma.if.usp.br/~konoplya/. This 
procedure (without any modifications) can be easily applied for calculations of characteris- 
tics of the Hawking radiation of various black holes. That is much easier than the shooting 
method which requires the analysis of asymptotical behavior of the wave equation for each 
black hole. 

One should note that for small to the accuracy of the WKB formula is small because of 
the large distance between the turning points in that case. For higher accuracy for small 
uj is achieved by the other formula, 



T = exp - / dn^V{n)-uj^ , (6.4) 



where the integration is performed between the two turning points V{r^) = uP' . Though 
the sixth-order WKB formula provides an excellent accuracy for the gravitational pertur- 
bations, for which £ > 2, even for small w. The matching of both WKB formulas, (|6.4| ) for 
small a; and ( |6.2| ) for moderate cj, provides usually very good estimation of the Hawking 
effect. 
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